arXiv:1502.03298vl [cond-mat.supr-con] 11 Feb 2015 


The Effect of Pointlike Impurities on Charge Density Waves in Cuprate 

Superconductors 

W. A. Atkinson^ and A. P. Kampf^ 

^Department of Physics and Astronomy, Trent University, Peterborough Ontario, Canada, K9J 7B8 
^ Theoretical Physics III, Center for Electronic Correlations and Magnetism, 

Institute of Physics, University of Augsburg, 86135 Augsburg, Germany 
(Dated: February 12, 2015) 

Many cuprate superconductors possess an unusual charge-ordered phase that is characterized 
by an approximate d„, 2 _y 2 intra-unit cell form factor and a hnite modulation wavevector q*. We 
study the effects impurities on this charge ordered phase via a single-band model in which bond 
order is the analogue of charge order in the cuprates. Impurities are assumed to be pointlike and 
are treated within the self-consistent t-matrix approximation (SCTMA). We show that suppression 
of bond order by impurities occurs through the local disruption of the d„, 2 _y 2 form factor near 
individual impurities. Unlike d-wa,ve superconductors, where the sensitivity of Tc to impurities can 
be traced to a vanishing average of the d„, 2 _y 2 order parameter over the Fermi surface, the response 
of bond order to impurities is dictated by a few Fermi surface “hotspots”. The bond order transition 
temperature Tbo thus follows a different universal dependence on impurity concentration Ui than 
does the superconducting Tc. In particular, Tbo decreases more rapidly than Tc with increasing m 
when there is a nonzero Fermi surface curvature at the hotspots. Based on experimental evidence 
that the pseudogap is insensitive to Zn doping, we conclude that a direct connection between charge 
order and the pseudogap is unlikely. Furthermore, the enhancement of stripe correlations in the 
La-based cuprates by Zn doping is evidence that this charge order is also distinct from stripes. 


I. INTRODUCTION 

Hole-doped cuprate superconductors have a pro¬ 
nounced “pseudogap” phase, which extends across a large 
fraction of the phase diagram. The physical origins of 
the pseudogap are unsettled, and the recent discovery 
of charge orderii^ in the pseudogap phase of a vari¬ 
ety of cupratesll*l^ has led to questions about a pos¬ 
sible relationship between the two. However, the con¬ 
nection is not straightforward: while the onset tem¬ 
perature for charge order T^o coincides with the tem¬ 
perature T* at which the pseudogap opens in single¬ 
layer Bi2Sr2_a;Laa;Cu06-i-5,^Tco is s ubstantia lly smaller 
than T* in other hole-doped cupratesj ^ l ^ l ^^b7 | j jg 
stantially higher than T* in the electron-doped cuprate 
Nd2_a:Cea;Cu04.1i^ It has been argued that some com¬ 
bination of charge, superconducting, and current fluctu¬ 
ations may persist up to T* and could be respon sible 
for the pseudogap in the hole-doped cupratcs.^^Hal Con¬ 
versely, some experiments appear to indic ate that charge 
order is distinct from the pseudogajP^IlIl and important 
aspects of the charge-ordered phase can be explained nat¬ 
urally und er th e assumption that it grows out of the 
pseudogap .iMUl 

The charge order has two distinguishing characteris¬ 
tics. The first is that it appears to have an approximate 
“nematic” or form factor, which is most easily 

understood as a transfer of charge between oxyge n sites 
along the x and y axes in the Cu02 planes.^^S®®! The 
strongest evidence for this comes f rom tunneling exper¬ 
iments in Bi-based cuprates, b l ^^ l ^^ l and it is further sup¬ 
ported by recent x-ray experiments.^^HnSl jg g^jgQ j-iote- 

worthy that form factor is widely predicted in 

calculat ions 


The second characteristic is that the amplitude of the 
interorbital charge transfer is modulated, with wavevec- 
tors q* = (q*,0) a nd q* = (0,q*) oriented along the Cu- 
O bond directions! ^ * ^ * ^ * ^^ ! The or ientation of q* has been 
hard to understand theoretically,f2lIlIl2lIlZIlll and in gen¬ 
eral calculations strongly prefer q-vectors oriented along 
the Brillouin zone diagonals. 

The charge order may thus be thought of qualita¬ 
tively as a “(i2,2_y2 charge density wave”. This charge 
density wave (CDW) appears to be qualitatively dif¬ 
ferent from the stripe order that has been widely ob¬ 
served in the La-base d cuprates La2-a;Baa:Cu04 and 
La2-a;Sra;Cu04.^ ^ * ^^ * ^^ l Perhaps the most compelling dis¬ 
tinction is that, whereas stripes in the La-cuprates have 
static or quasistatic spin and charge modulations whose 
periods are locked together, there is no apparent cor¬ 
relation between spin and charge degrees of freedom in 
YBa2 Cu3 Oe-i-x 

In this work, we address the question of how a dx 2 -y 2 
CDW responds to strong-scattering pointlike impuri¬ 
ties. Furthermore, because charge order is known to 
coexist with superconductivity at low temperatures,!^ 
we explore the effects of impurities on a mixed 
superconducting-CDW phase. We adopt a simplified 
one-band model in which the analogue of charge order is 
an anisotropic renormalization of the electron hopping, 
known as bond order. Bond order and superconductiv¬ 
ity are driven by a combination of spin exchange and 
Coulomb interactions between nearest-neighbor lattice 
sites.!^ 

While the structure of t he C DW will be affected by 
any preexisting pseudogap) ^^ * ^^ ! we avoid complications 
associated with modeling the pseudogap and assume in¬ 
stead that the CDW grows out of the full Fermi surface. 























2 


The impurity physics described in this work is sufficiently 
general, however, that it should equally apply in the pseu¬ 
dogap phase. 

We distinguish here between two separat e issu es. First, 
it has been pointed out by several authord^Ulz] 
directional charge density waves break both a C/(l) sym¬ 
metry associated with the location of the CDW and a 
Z2 symmetry associated with its orientation. Disorder 
couples linearly to the CDW, and immediately restores 
the C/(l) symmetry (in the disorder average), leaving an 
“electron nematic” phase that breaks rotational but not 
translational symmetry. In the nematic phase, the model 
then maps onto the random field Ising model, which in 
three dimensions has a critical disorder strength above 
which long range rotational order is destroyed. We note, 
however, that while long range CDW or nematic order is 
destroyed, local CDW order persists on a length scale set 
by the disorder potential.!^ 

The second issue concerns the suppression of the ampli¬ 
tude of the charge order by impurities. X-ray scattering 
experiments on YBa2Cu306.6 observe a rapid reduction 
of the charge order by Zn impurities.!^ Naively, this is 
expected given the dx2_y2 form factor of the charge or¬ 
dered state: any order parameter whose average over the 
Fermi surface is zero should be rapidly suppressed by 
isotropic scattering. This mechanism is responsible, for 
example, for the well-known breakdown of Anderson’s 

theorem in d-wave superconductorsgg pointed 

out previously by Ho and Schofield,!^ a mathematically 
identical theory describes the suppression of the second 
order d^^-y^ Pomeranchuk Fermi surface instability. 

Here, we show that the d^^-y^ symmetry of the charge 
order is of marginal importance; rather, there is a rapid 
suppression of charge order by impurities that can be at¬ 
tributed to the “hotspot” structure of the charge ordered 
phase. In particular, the rate at which charge order is 
suppressed depends sensitively on the Fermi surface cur¬ 
vature near the hotspots. Consequently, the suppression 
of bond order by impurities follows a different universal 
relationship than d-wave superconductors. 

This is of direct relevance to the cuprates, where 
Zn substitutes isovalently for Cn an d acts as a strong¬ 
scattering pointlike impurity! ^^ * ^^ ! "*^^ ! Zn-doping has, in 
past, been used as an important local probe of the 
superconducting^ and pseudogap states,! 2 ^M^. Our work 
leads us to three main conclusions: first, we find that 
charge order is more rapidly suppressed than d-wave su¬ 
perconductivity; second, the insensitivity of the pseudo¬ 
gap to Zn doping makes charge order an unlikely cause 
of the pseudogap; third, the insensitivity of stripes in 
the La-based cuprates to Zn doping supports that the 
stripe physics is inherently different from charge order in 
YBa 2 Cu 306 +x and the Bi-based cuprates. 

We begin in Sec. miby introducing the model, and pro¬ 
vide an introductory discussion of how bond order re¬ 
sponds to pointlike impurities. Disorder-averaged equa¬ 
tions for the bond order parameter are then derived for 
the case of a dilute concentration of strong scattering 


impurities, treated within a self-consistent t-matrix ap¬ 
proximation (SCTMA). 

In Sec. |HI[ we consider temperatures near the bond or¬ 
der transition temperature Tbo where the equations can 
be linearized. These equations are the same for uni¬ 
directional and bi-directional (checkerboard) order be¬ 
cause the different Fourier components of the bond order 
decouple near Tbo, and their simplicity allows one to de¬ 
rive approximate analytic expressions for the dependence 
of Tbo on the impurity concentration rii. We find that 
Tbo is suppressed more quickly by disorder than is the 
transition temperature Tc for d-wave superconductivity. 

In Sec. IV A we obtain numerical mean-field solutions 
for T < Tbo for situations in which the bond order is 
commensurate with a periodicity of m unit cells. X-ray 
experiments on YBa2Cu306+a; generally find two sets of 
peaks, rotated by 90 ° relative to each other. The relative 
intensity of the peaks is strongly doping dependent!^ and 
in YBa2Cu3 06.54, where charge order is strongest, the 
peak along the a-axis is almost undetectably weak.l^ This 
suggests that the two Fourier components of the charge 
order are not strongly coupled. Furthermore, analyses 
of STM experiments in Bi-based cuprates find nanoscale 
doma ins of unidirectional order rather than true biaxial 
orderWe therefore focus on unidirectional order, 
although an extension to multiple Fourier components is 
conceptually straightforward. Our calculations lead to a 
coupled set of equations for the impurity scattering rate 
(or, more specifically, the self-energy) and the bond order 
parameters. Most of the technical details are relegated 
to the appendices. 

Experimentally, the charge ordering temperature T^o 
is greater than Tc, and we therefore examine the on¬ 
set of superconductivity in the presence of bond order 
in Sec. |IVB| We derive superconducting Tc equations 
within the SCTMA, which are then solved in conjunc¬ 
tion with the self-consistent equations for the bond order 
parameter. We find that, in the mixed phase, impurities 
suppress Tbo more rapidly than Tc, and one may therefore 
obtain Tc > Tbo as the impurity concentration increases. 

Finally, in Sec. |Vj we discuss our results in the context 
of Zn-doping experiments in the pseudogap phase, and in 
the stripe phase of La2-aBaa;Cu04. These experiments 
show that the pseudogap and stripe phases respond dif¬ 
ferently to Zn impurities than the charge ordered phase 
seen in YBa2Cu306+a;, and therefore likely have a differ¬ 
ent origin. 


II. MEAN-FIELD EQUATION FOR BOND 
ORDER 

Following Ref. EQI we adopt a tight-binding model on 
a square, two-dimensional lattice, representing a single 
Cu 02 plane. The noninteracting part of the Hamiltonian 
is 

Ho = CkC^cCka (1) 

k,Q; 
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with dispersion ek = to ~ 2ti{coskx + cosky) + 
At 2 cos kx cosky. We take ti = 1, which sets the energy 
scale for the calculation; in the cuprates the bandwidth 
8 ti is of order a few electron volts. The interacting part 
of the Hamiltonian contains both nearest-neighbor ex¬ 
change and Coulomb interactions, 

{i,j)a,...,S a (ij) 

( 2 ) 

where (z, j) refers to nearest-neighbor lattice sites i and 
j and are Pauli matrices. We perform a mean-field 
decomposition of this interaction in the exchange channel 
to obtain 

Hlho = ( 3 ) 


where is assumed to be independent of the spin 

a. This term drives the bond-ordering instability, and we 
therefore define an effective interaction for bond order. 



jso={^-^+vy ( 4 ) 


A similar decomposition in the particle-particle channel 
leads to a mean-field superconducting contribution, 


H 


Isc 


— J?,l 


H + h-c. 

(*b> 


( 5 ) 


with 


= (x - 1^) • (6) 

While the spin-exchange interaction is attractive in 
both the bond-order and superconducting channels, the 
Coulomb interaction enhances bond order and suppresses 
superconductivity. This was invoked previously as an ex¬ 
planation for why Tbo is greater than the superconduct¬ 
ing transition temperature 

We denote the exchange self-energy along the nearest- 
neighbor bond i-j by 

Pji = ~Jho{c^ja^i(y) ■ C^) 

For illustrative purposes, we solve this self-consistently 
in real space: Pji is calculated for each bond by diag¬ 
onalizing the mean-field Hamiltonian Hq + i?ibo on an 
Lx L lattice with periodic boundary conditions; the new 
Pji are used to update i?iboi and the process is iterated 
until self-consistency is achieved. To obtain a solution, it 
is necessary to tune the band parameters such that the 
system size is an integer multiple of the CDW period. 

Typical results for a clean lattice are shown in 
Fig. s a). To highlight the bond ordering, we have 
subracted off the uniform exchange self-energy Pq of 


FIG. 1: Results of self-consistent real-space bond-order cal¬ 
culations on a 48 X 48 lattice, (a) Bond self-energies Ptj be¬ 
tween nearest-neighbor lattice sites in the clean (impurity- 
free) limit. The constant reference self-energy Pq of the ho¬ 
mogeneous phase is subtracted, as described in the text, to 
highlight spatial fluctuations. Blue (red) line colors indicate 
that the self-energy is enhanced (reduced) relative to Pq, and 
linewidths indicate the size of the enhancement (reduction), 
(b) k-dependent spectral function at the Fermi energy for the 
clean limit, (c) Bond self-energies for a system containing a 
single impurity, indicated by a black circle at the center of 
the figure. Here, Po has not been subtracted when determin¬ 
ing the linewidths; however, the color scheme is the same as 
in (a), (c) Bond self-energies for a concentration rii = 0.04 
of impurities. Impurity locations are indicated by black cir¬ 
cles. Model parameters are to = 0.85, ti = 1.0, t 2 = —0.5, 
Jbo = 2.40, Vi = 10.0, and T = 0.25. 


the homogeneous phase, which is obtained by requir¬ 
ing all bonds to be equivalent in the self-consistent cal¬ 
culation. The amplitudes and signs of the shifts in 
the bond self-energy, relative to Pqi shown by the 
thicknesses and colors of the lines connecting nearest- 
neighbor sites. For the parameters chosen, the modula¬ 
tion amplitude is about 10% of Pq. The phase shown in 
Fig. I^a) has Fourier components of equal magitude at 
q = ±(27r/4)(l, 1) and q = ±(27r/4)(l, —1), and can be 
thought of as a dx^-y^ form factor whose amplitude has 
a period-4 checkerboard modulation. 

The corresponding spectral function A(k,£i?) at the 
Fermi energy is shown in Fig. 0b)- The bond-order q- 
vectors connect segments of Fermi surface near (±7r, 0) 
and (0, ±7r), and which points a spectral gap opens. The 
residual Fermi surface has a strong intensity along arcs 
centered on the Brillouin zone diagonals. Faint residual 
Fermi surface segments can also be seen along the Bril- 
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louin zone boundaries. 

Figure [^c) shows the effects of adding a single strong¬ 
scattering impurity to the lattice, modeled as a potential 
shift of Vi = 10.0 at position (25,25). In this panel, the 
linewidths indicate the total bond self-energy, including 
Pq; however, the color scale still indicates the shift of 
Pij relative to Pq, as in Fig. [ija). Unsurprisingly, 
is reduced almost to zero along bonds connecting to the 
impurity site. The amplitude of the bond self-energy re¬ 
covers within a lattice spacing of the impurity; however, 
the modulation pattern is disrupted over a longer length 
scale. Thus, suppression of bond order does not simply 
imply that the bond self-energies vanish, but rather that 
the form factor is disrupted near impurities. This is ev¬ 
ident in Fig. 0d) where the structure shown in (a) is 
entirely disrupted by a 4% impurity concentration. This 
scenario is to be contrasted with that presented in, for 
example. Ref. where smooth impurity potentials dis¬ 
rupt long range order but leave the form factor locally 
intact. 

To proceed further, we study the disorder-averaged 
mean-field equations for bond order, and subsequently 
superconductivity, for a dilute concentration of strong 
scattering impurities. Because disorder-averaging re¬ 
stores translational symmetry, it is useful to Fourier 
transform Eq. 0 to k-space, 

^Ibo — ^k(q)c],+q^Cka (8) 

k,q,tT 

where 

^k(q) = (9) 

i,5 

= E 2 [cos(A:^ - K) + cos{ky - k'y)] 

k' 

^ ('"k'o-'^k'-l-qCT)- (10) 


with 


77k = V2cos{kj; + y), 

(12a) 

vl = V2cos{ky + ^), 

(12b) 

tjI = y2sin(fc:r + y), 

(12c) 

-f 

II 

(12d) 

reduces to 


Tk(q)=E<^“(q)’ 

(13) 


with 


^’“(q) = -^E'?k'(Ck'<rCk'-rqa). 
k^ 

Jho 


(14) 


= -^E^k'G(k'-f q,k';T = 0 ), (15) 


where G'(ki,k 2 ;r) is the Green’s function at imaginary 
times r in the presence of bond order. This equation has 
an even solution 


^’k(q) = -pi(q)?7k + -P^(q)^L 

and an odd solution 

i^k(q) = U'(q)r?^ + U"(q)^^ 


(16) 


(17) 


The even solution is the leading instability in all calcula¬ 
tions reported here. 

Equation (15) requires an explicit expression for the 
Green’s function, and we consider two cases where closed 
expressions are possible: (i) temperatures near Tbo where 
a linearized Green’s function can be obtained and (ii) the 
case of period-TO commensurate order. 


III. LINEARIZED EQUATIONS FOR Tbo 


In this definition of Pk(q), k and k-l-q are initial and final 
momentum labels for electrons that are scattered by the 
bond order; a common alternative is to take these initial 
and final points to be at k —q/2 and k-l-q/2. This latter 
choice is less convenient for systems with commensurate 
bond order, which we discuss below. 


Equation (10) is the basic self-consistent equation for 
Pk(q)- It is invariant under k ——k — q when Pk(q) 
is real, so that Pk(q) should be even or odd under this 
transformation. We then expand Eq. (101 in a set of basis 
functions that are even or odd under k —>■ — k — q, via 


2 [cos(fca; 


fc(,) -I- cos{ky 


4 

= E^“(k),7“(ko, 


q;=1 

( 11 ) 


Near the bond ordering transition, Pk(q) is small, so 
that iJ ibo c an be treated perturbatively. We show in Ap- 
pendix [AT] that to linear order in Pk(q), Eq. (15) reduces 
to a matrix equation for the components P“(q), 


p“(q) = jboE^“/5(q)^^(q)- (18) 

p 


Equation ( |18| ) has a solution at the ordering wavevector 
q* when the largest eigenvalue of the matrix F(q*) is 
equal to 1/Jbo- For a given Jbo, we search for the tem¬ 
perature Tbo at which bond order emerges. We note that 
the different q-vectors are decoupled at Tbo, with each 
independently satisfying an equation of the form (181. 
Gonsequently, the dependence of Tbo on impurity con¬ 
centration is the same for uni-directional charge order as 
it is for bi-directional (checkerboard) order. 
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FIG. 2: Comparison of impurity effects on period-4 axial bond 
order and d-wave superconductivity. The bond ordering tem¬ 
perature Tbo and superconducting transition temperature Tc 
are shown as a function of the impurity concentration for im¬ 
purity potential Vi = 10. Results are shown for full numerical 
calculations (solid curves), and for the analytical equations 
( |23| l and ( |24[ | (dashed curves). We have taken k = 1.0 in 
Eq. (23l. To determine 7 for the analytical equations, a lin¬ 
ear ht was made to —Im for n = 0 , as a function of rn: the 
resulting formula is given in the figure. Inset: Fermi surface 
and q*. Band parameters are to = 0.85 and t 2 = 0.5. The 
interaction strengths Jbo = 2.74 and Jsc = 1.44 are chosen to 
obtain = r“ = 0.15. 


A plot of Tbo versus impurity concentration rii is shown 
in Fig. for period-4 axial order, with q* = (27r/4,0), 
for strong scattering impurities. To obtain this figure, we 
have tuned the band parameter to, which controls the fill¬ 
ing, such that q* connects parallel segments of Fermi sur¬ 
face (see Fig. inset). For comparison, the depen denc e 
of To on rij, calculated with the SCTMA (Appendix |B 1[ ), 
is shown for a d-wave superconductor. We see that both 
bond order and superconductivity are suppressed by dis¬ 
order, but that Tbo is suppressed more rapidly than T^- 
To understand the suppression of bond order by im¬ 
purities, we analyze the eq uatio ns governing Tbo. The 
kernel TQ,^(q) is (Appendix |A 1[ ) 

F fgV ^ dk(dk +‘^n’^) 

N ^ - Ck+q - s„)(ia;„ - ek - S„) ’ 

(19) 

where a;„ = (2n-|- l) 7 rT are Matsubara frequencies, is 
the impurity self-energy at uin to zeroth order in Tk(q), 
and the linear-order impurity self-energy S)* has been 
factored into components, 

S^ = ^5yT^(q). (20) 

/3 


For temperatures near Tbo, the SCTMA gives the self- 
consistent equation for the zeroth order impurity self- 


energy [Eq. (|A4|], 
Yin — 


i-Yiy 

N 


^k ^r; 

k 


( 21 ) 


where rii is the impurity concentration and V) is the im¬ 
purity potential. The real part of acts as a chemical 
potential shift due to doping by the impurities, and the 
imaginary part is the negative of the scattering rate 7 „. 
Because 7 „ has the same sign as uJn, it behaves qualita¬ 
tively like a temperature increase: it can be absorbed into 
a renormalized Matsubara frequency (!>„ = u!n+"/n whose 
magnitudes are larger than the unrenormalized frequen¬ 
cies ojn- The effect of is therefore to reduce bond 
order and suppress Tbo. 

The physics of is quite different from that of I]„. 
We find numerically that when is omitted from the 
self-consistent calculations, Tbo is reduced. This is simi¬ 
lar to the situation in superconductors, which we review 
in Appendix |B ![ where an analogous “anomalous” self¬ 
energy appears in the equations for T^. In conventional 
isotropic s-wave superconductors, the Tc enhancement by 
the anomalous self-energy cancels the reduction of Tc by 
En [cf. Eq. (B20)], consistent with Anderson’s statement 
that Tc is unaffected by disorder. E3 The response of Tc to 
impurities is closely tied to the symmetry of the super¬ 
conducting order parameter: in d-wave superconductors, 
the anomalous self-energy vanish es [cf. Eq. (BI 6 I] and T^ 
is strongly re duce d by impurities . 1^21111 

From Eq. (A9), the expression for E)* is proportional 


to a weighted average of Tk(q) over the Brillouin zone: 


E3oc 


^k(q) 


k ^k ^k+q Ey,) 


( 22 ) 


The sum in Eq. ( |22[ ) is weighted towards those points, the 
so-called “hotspots”, for which k and k -|- q both lie on 
the Eermi surface. This has two consequences: first, the 


k-sum in Eq. (22) does not vanish, even when Tk(q) has 
a nominally d-symmetric form factor Tk(q) oc ??k ~ ^k> 
second, E)* nonetheless tends to be small because of the 
limited region of k-space that contributes to the sum. 
Indeed, the omission of E)J from the self-consistent cal¬ 
culations only changes Tbo by a few percent. We conclude 
that the sensitivity of Tbo to impurities is not tied to the 
symmetry of the form factor, but is a consequence of the 
central role of hotspots in the Tbo calculation. 

We can integrate Eq. (19) analytically under a few sim¬ 


plifying assumptions. We expand the electronic disper¬ 
sion around the Fermi surface hotspots, and ignore the 
energy dependence of E„, letting E„ —>• —i 7 Sgn(w„). We 
obtain (see Appendix | A 2\ 


Tbo 


dkyKe 




,( 23 ) 
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where k is the Fermi surface curvature at the hotspots, 7 
is the scattering rate, is the bond ordering tempera¬ 
ture in the clean limit, and 'ijj{x) is the digamma function. 
Equation ( |23[ ) obtains a form similar to the usual result 
for the transition temperature of a d-wave superconduc¬ 
tor, namelj®^ 


'T'Q 



(24) 


when the Fermi surface curvature is At = 0. 

These analytical expressions are shown in Fig. To 
make the comparison quantitative, we have set 7 = 
—Im where ujq is the lowest positive Matsubara 

frequency. The curvature is 


K = 




dk^ 


k=(27r— ^ ,7r) 


which is typically a number of order 1. It is apparent 
in Fig.j^that the analytical expressions overestimate the 
transition temperatures somewhat, but that they capture 
the reduction of Tbo relative to Tc. We note that the sen¬ 
sitivity of Tbo to disorder is in addition to the reduction 
of TiCq due to K in the clean limit; indeed, in Fig. |^we 
had to use an inflated value of Jbo = 2.74 relative to the 
pairing interaction Jgc = 1.44 to obtain T^^ = T°. 


IV. COMMENSURATE BOND ORDER 

A. Pure bond order 


When the bond order parameter is not small, we 
can proceed by assuming that the wavevector is com¬ 
mensurate, with q* = (27r/TO)(l,0) for axial order and 
q* = 27r/TO(l,l) for diagonal order, where m is an in¬ 
teger. These describe uni-directional phases, and the 
extension to bi-directional order is straightforward. For 
clarity, we describe only the case of uni-directional order. 

When the bond modulation has a period of m unit 
cells, the Brillouin zone is correspondingly reduced by a 
factor of m along one direction. This is illustrated in 
Fig-i The mean-field Hamiltonian = Ho + Hioo 
can be written in matrix notation as 

= E E^bHk(q*)«'k. (25) 

kGBZ' cr 


where BZ' is the reduced Brillouin zone, and Tk is a col¬ 
umn vector of length m containing annihilation operators 
with momenta connected by integer multiples of q*: 


«'ka 


Cki cr 
Ck 2 CT 


(26) 



FIG. 3: Construction of the reduced Brillouin zone. In this 
example we consider period-4 axial order (m = 4), with mod¬ 
ulation wavevector q* = (27r/4)(l,0). The original Fermi 
surface (solid black curve) is shown in the full Brillouin zone, 
which extends from 0 to 2n along and ky. The reduced 
Brillouin zones have width |q*| aloud kx, and extend from 0 
to 27r along ky. The different reduced Brillouin zones are la¬ 
beled j = 1,2,... ,m, and the principal reduced Brillouin zone 
(shaded region) has j = 1. We show the zone-folded Fermi 
surface in the principal reduced zone. Any point ka belonging 
to the ath reduced zone can be written = k 4- (a — l)q*. 
Two such points, k and k 2 are illustrated in the figure. 


with 


ka = k-h (a-l)q*, ae[l,m]. (27) 

In this notation, ka belongs to the ath reduced Brillouin 
zone. The m x m matrix Hk(q*) has nonzero elements 


{ ek„, a = b 

Fk 6 (q*), a= mod ( 6 , m)-hi . 
Pk„(q*)*, b= mod (a, m)-h 1 

(28) 

Then, the matrix Green’s function (including the im¬ 
purity self-energy matrix S„) is G(k, fa;„) = — 

Hk(q*)-Sa]-1 with matrix elements 

Ga6(k; u}) = G(ka, kf,; w). (29) 


Substituting this into Eq. (15), the equations for the bond 
order follow: 


T“(q*) 


•^bo 

ruN' 


kSBZ' a=l 


0-) (30) 


where it is understood that a -I- 1 = mod (a, m) -I- 1 
and N' = N/m is the number of k-points in the reduced 
Brillouin zone. 

Without disorder, we can evaluate the Green’s func¬ 
tion from the eigenvectors and eigenvalues of Hk(q*) to 
obtain 

J m 

^“(q*) = -;^ E E C^“+i/(k)^a7(k)/(i?ik) 

kGBZ' a,e=l 

(31) 


Ck^(7 
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FIG. 4: Temperature dependence of the order parameter for 
the q = (0, 0) nematic instability. Results are shown for con¬ 
centrations Ui of strong scattering pointlike impnrities with 
potential Vi = 10.0. Inset. Comparison of the nematic Fermi 
surface (red solid) and the bare Fermi surface (black dashed). 
Model parameters are to = 1.55, = 0.4, and Jbo = 1.5. 


where S(k) is the matrix of eigenvectors of Hk(q*), and 
Eii^ are the corresponding eigenvalues. More generally, 
once disorder is included, we have 


G(k;r = 0") =r 




[tw„l - Hk(q*) - S, 


1-1 


(32) 

Substitution of Eq. (32) into Eq. (301 generates the self- 
consistent equation that must be solved for P“(q*). The 
prescription for obtaining within the SCTMA is de¬ 
scribed in Appendix |A 3| 

As a point of reference, we first revisit the case of ne¬ 
matic order (ie. the Pomeranchuk instability) which was 
previously studied by Ho and Schofield for a Gaussian 
distributed disorder potential.!^ The Pomeranchuk tran¬ 
sition is a q = 0 instability (so m = 1), and to ob¬ 
tain it one must tune the Fermi surface so that it passes 
near the Brillouin zone boundaries at (± 7 r, 0) and (0, ± 7 r). 
Here, we take the next-nearest neighbor hopping ampli¬ 
tude t 2 = 0.4, and adjust the filling to obtain the Fermi 
surface shown by the dashed curve in the inset of Fig. 

The leading q = 0 instability has a pure dj; 2 _y 2 (or 
nematic) symmetry, with 


Pk{0) = P[cos{kx) - cos{ky)], (33) 

where P = ^pi(q *)2 _|_ p 2 ^q*^ 2 ^ resulting Fermi 

surface in the bond ordered phase is shown by the solid 
curve in the inset to Fig.|^ the Fermi surface distortion 
has a clear dx 2 _y 2 symmetry, with points near ( 0 ,± 7 r) 
pushed in and points near (±7r, 0) pushed away from the 
Brillouin zone center. 

Impurities suppress P, and thereby this distortion, as 
shown in the main panel of Fig. and the nematic 
phase is ultimately destroyed near rii « 0.10. Ho and 



FIG. 5: Phase diagram for period-5 diagonal order, with 
q* = (27r/5)(l, 1). (a) Original (dashed red) and recon¬ 

structed (solid black) Fermi surfaces. The arrow labelled q* 
connects Fermi surface hotspots that are gapped by bond 
order, (b) Onset temperature as a function of interaction 
strength for clean {rn = 0.00) and disordered (n^ = 0.04 
and Ui = 0.08) systems with strong scattering impurities 
{Vi = 10.0). (c) Magnitude of the order parameter as a func¬ 
tion of temperature for Jbo = 1.68. Model parameters are 
to = 1.3 and t 2 = 0.5. 


SchofielcP^ noted previously that disorder can change the 
order of the transition from second to first in cases where 
the Fermi surface does not pass exactly through (± 7 r, 0) 
and (0, ± 7 r). This same crossover can be seen in Fig. 
at rii « 0.04. In cases where the nematic transition is 
second order, Tbo satisfies the same dependence on the 
impurity scattering rate 7 as d-wave superconductivity,!^ 
namely Eq. (24). 

Next, we examine the case of diagonal order which, as 
discussed in the Introduction, was the leading instability 
in a large number of earlier calculations. We choose q* = 
(27r/m)(l,l) with m = 5 to give q * that i s similar in 
magnitude to what was found earlier To obtain a 
solution, it is necessary to tune the band parameters so 
that q* connects antiparallel hotspot sections of Fermi 
surface, as shown in Fig. [^a). Near the hotspots, bond 
order gaps the Fermi surface and thereby reconstructs it 
as shown in Fig. [^a). 



















We remarked earlier that the self-consistent equation 
for Pk(q), Eq. (10), is invariant under k — )■ —k — q. For 


diagonal order, Eq. (10) is also invariant under ky. 


Based on these two symmetries, we expect solutions for 
Pk(q*) to have the form 


^k(q*) = P 


* 

cos{k^ + Y)^ cos{ky + -^) 


(34) 


In our calculations, the solution with the negative sign 
is always preferred. While this solution superficially re¬ 
sembles the d-symmetric order parameter found at q = 0, 
Pk(q*) does not have even a qualitative interpretation as 
a dx 2 -y 2 distortion of the Fermi surface. Indeed, because 
of the Brillouin zone folding associated with the finite- 
q modulation, the reconstructed Fermi surface shown in 
Fig. §a) is quite complicated, with no resemblance to 
that in Fig. 

We show the dependence of Tbo on Jbo for different im¬ 
purity concentrations in FigjMb), and the T-dependence 
of P for different in Fig. ^c). Similar to the nematic 
transition, impurities reduce P; here, however, the ne¬ 
matic transition remains second order as the impurity 
concentration grows. As in Sec. |III[ the different q com¬ 
ponents of the order parameter decouple near Tbo, and 
Tbo is the same whether the order is uni-directional or 
bi-directional (checkerboard). 

Finally, we consider axial order with q* = 
(27r/TO)(l,0), as shown in Fig. We take m = 3, 
which gives q* close to that seen experimentally in 
YBa2Cu306+a;- Again, it is necessary to tune the band 
parameters such that q* connects antiparallel portions of 
the Fermi surface [Fig. [^a)]. To enhance the suscepti¬ 
bility towards axial order, we have taken t 2 = 0.7, which 
reduces the curvature near the Fermi surface hotspots. 
(The connection between curvature and Tbo is discussed, 
e.g. in Ref. |34l) Nonetheless, a rather large Jbo = 3.8 
is required to obtain a clean-limit transition temperature 
T°q that is the same as in Fig. |^for the nematic insta¬ 
bility. 

In the axial case, the self-consistent equation for 
Tk(q*) is invariant under k^ ^ —kx — q* and ky —> —ky. 
This implies that the order parameter has the form 


E’k(q*) = V2T 


cos(a) cos{kx fo fo sin(a) cos{ky) 


(35) 

where tana = P^(q*)/T^(q*) is a nonuniversal constant. 

Figure [^b) shows Tbo as a function of Jbo for two 
different values of n^, along with tana. As before, these 
results hold for both uni-directional and bi-directional 
order. From the plot of tana, we see that the magnitude 
of P^(q*) is 60-70% of the magnitude of P^(q*), that 
P^(q*) has the opposite sign of P^(q*), and that disorder 
changes this admixture. 

It is notable that the bond order in the axial and diag¬ 
onal cases is more rapidly suppressed than in the nematic 
case, with the axial case the most sensitive to impurities. 
Equation (23) suggests that in the axial and diagonal 



X 




FIG. 6: Phase diagram for period-3 axial order, with q* = 
(27r/3)(l,0). (a) Original (dashed red) and reconstructed 

(solid black) Fermi surfaces. The arrow labelled q* indicates 
the bond ordering vector connecting Fermi surface hotspots, 
(b) Phase boundaries (solid symbols) for the onset of axial 
bond order for clean (n; = 0.00) and disordered {rii — 0.04) 
systems with strong scattering impurities {Vi = 10.0). The 
relative amplitude tana of the ky and components of the 
order parameter, defined by Eq. ( |35| ), is also shown (open 
symbols), (c) Temperature dependence of the order parame¬ 
ter P at different impurity concentrations for Jbo = 3.8. The 
dependence of a on T is very weak. Model parameters are 
to = 0.4, t2 = 0.7. 


cases, Tbo depends on both the Fermi surface curvature 
and scattering rate. The nematic transition, on the other 
hand, appro xim ately satisfies an equation of the same 
form as Eq. (24) for J-wave superconductivity,!^ and at 
this level of approximation depends only on the scat¬ 
tering rate; nematic order is thus expected to be more 
robust against impurities than finite-q bond order, con¬ 
sistent with the numerical results shown in Figs. Hi and 
[^ Furthermore, comparing the diagonal and axial cases, 
we note that the Fermi surface curvature in the diagonal 
case {k « 0.3) is approximately half that for the axial 
case {k « 0.6), which is consistent with the more rapid 
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suppression of Tbo in Fig. 1^ than in Fig. 

Importantly, the scattering rate 7 also depends on 
band structure. In the strong-scattering limit (Vi —>■ 00 ), 


7 = 


TTUi 


(36) 


where Nq is the density of states at the Fermi energy. 
The scattering rate is thus smallest for the nematic order 
in Fig. [^because the Fermi surface passes near van Hove 
singularities at (± 7 r, 0) and ( 0 ,± 7 r). We find that for a 
fixed rii the scattering rate for the axial case is roughly 
twice that for the nematic case, and slightly less than 
twice that for the diagonal case. These differences in 7 
are consistent with the different sensitivities to impurities 
shown in Figs. 00 and 1^ In summary, the sensitivity of 
bond order to impurities depends on the band structure, 
both directly through the Fermi surface curvature and 
indirectly through the scattering rate. In cuprates, we 
can thus expect that the sensitivity of charge order to 
impurities will be doping-dependent. 


B. Tc equations for superconductivity in the bond 
ordered phase 


To explore the onset of superconductivity in the bond 
ordered phase, we consider linearized equations for the 
pairing instability in the presence of period-m commensu¬ 
rate bond order. These will give both the superconduct¬ 
ing transition temperature Tc, and the k- and q-structure 
of the order parameter Ak(q) near Tc- The mean- 
field pairing contribution to the Hamiltonian, Eq. (§ is 
Fourier transformed to obtain 

H\sc — X! [^k(q)*c_k+qiCkt + h.c.], (37) 

k.q 

with Ak(q) = Ea 7 £A“(q) and 

A“(q) = -Jsc^ Xl^?k'(c-k'+q;Ck't). (38) 

k' 

The basis functions defined by Eqs. ( [T^ , are the 
same as used to describe the bond order. If the bond 
order has wavevector q*, then the pair order parameter 
A“(q) must necessaril y hav e Fourier components q = 
0, q*, 2q*,... (m — l)q* .ISHISI Defining 


= 


‘^kit’ ''k2T’ ■ ■ ■ ‘"bf '^-kit’ '^-k2ti 


(39) 

where = k -|- (a — l)q*, the mean field Hamiltonian 
containing both superconductivity and bond order is 


Hm{= ^ 


kGBZ' 


Hk(q*) 


A 


-H^k(-q*) 


^k, (40) 


where Hk(q*) is defined in Eq. (28), and the mx m off- 
diagonal block has matrix elements 


[A^U = Y.^a-bV^.- ( 41 ) 

Ct 


In this expression, A“_^ is shorthand for A“[(a — &)q*]. 
The diagonal elements [Ak]aa therefore correspond to 
pairs with zero center-of-mass momentum belonging to 
the ath reduced Brillouin zone. We note that because 
mq* is a reciprocal lattice vector we use terms like (m — 
l)q* and —q* interchangeably. 

The expectation value in Eq. ( [M| can be evaluated to 
linear order in the pair amplitude to obtain the eigenvalue 
equation for the elements of Ak (see Appendix |B 2[ ) 

= Jsc E E M^ia-byMc-d)Atd (42) 

C,d 

with 

Maa-fic = E X! ^kf^kf, 

keBZ' n = l 

X \iujji Hk(q ) 

x[za;„ + H^k(-q*) + S:]Ec7-a- (43) 


In this equation, it is understood that (! — c and i — a 
are evaluated modulo m. Furthermore, we have dropped 
the anom alous impurity self energy as discussed in 


Sec. 


HI 


Tin vanishes identically in pure d-wave super¬ 


conductors, and as we show below, superconductivity 
has predominantly d-wave symmetry in the bond-ordered 
phase. The neglect of leads us to underestimate 
Tc slightly; however, there are two relevant cases where 
Tn = 0 exactly: (i) rii = 0 , where the impurity self en¬ 
ergy vanishes, and (ii) cases in which impurities suppress 
Tbo such that Tc > Tbo, and the superconductivity is 
purely d-wave. 

The kernel Maa-,pc forms a 16 x 16 matrix M with rows 
and columns labeled by the composite indices {a, a) and 
(/3, c) respectively. The superconducting instability oc¬ 
curs when the largest eigenvalue of M is equal to 1/dsc- 
We show the dependence of Tc on impurity concentration 
in the axial bond-ordered phase in Fig. The bond or¬ 
der parameter Tk(q*) is calculated self-consistently, and 
Tbo is therefore also shown. Experimentally, charge order 
emerges at a higher temperature than superconductivity, 
although the ratio of Tc and Tco is doping dependent and 
decreases with increasing hole concentration.lISl We there¬ 
fore show two cases in Fig. 0 In the first, T^ = T^o/3, 
which is comparable to the smallest ratio of Tc to the 
charge-ordering temperature seen by x-ray experiments 
in YBa 2 Cu 3 06-i-a:-^ As Ui increases, Tbo decreases faster 
than Tc, although superconductivity is destroyed hrst. 
In the second case, T° = 2T^^/2,, which is slightly larger 
than the maximum ratio found in YBa 2 Cu 306 +a;. Here, 
there is a narrow window over which impurities destroy 
bond order, but superconductivity remains. Note that 
our calculations explicitly neglect the feedback of super¬ 
conductivity on bond order, and therefore overestimate 
Tbo in regimes where Tbo < Tc (although Tc is correctly 
given). 

The eigenvector corresponding to the largest eigen¬ 
value of M gives the k- and q-space structure of the 
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(c) X (d) X 

FIG. 7: Superconducting transition in the bond-ordered 
phase. Results are shown for period-4 axial bond order, (a) 
Transition temperatures for bond order and superconductiv¬ 
ity. Results are shown for Jbo = 2.74 and two values of the 
pair interaction Jsc- Model parameters are otherwise the same 
as in Fig. (b) Fourier components of the superconducting 
order parameter at Tc, relative to the d-wave component vq 
[cf. Eq. (|45[)], for Jsc = 1.31. (c) Real space pattern of the 
nonzero-q components of Pij at rii = 0, T — 0.10. The order 
parameters P^(q*) = 0.102 and P^(q*) = —0.082 are self- 
consistently determined, (d) Nonzero-q component of the su¬ 
perconducting order parameter Aij at rii — 0 at T = 0.10. 
The pattern is derived from the pairing kernel eigenvector 
given in Table |T] In (c) and (d), the line thickness and color 
indicate the magnitude and sign of the order parameters re¬ 
spectively. Band parameters are to = 0.85 and t 2 = 0.5. 


electron pairs near the superconducting transition. Given 
the eigenvector Vaa, we then have 


Ak(q) = A(T) ^ l5q,(a-l)q* 

a—1 


■ 4 

_Q !=1 


(44) 


where A(T) is the amplitude of the order parameter and 
the terms in the square brackets give the A:-space struc¬ 
ture of each Fourier component of Ai(;(q). This equation 
makes explicit that the pair wavefunction has contribu¬ 
tions at multiple center-of-mass momenta. 

To give a concrete example, we consider the order 
parameter for = 0 at temperatures slightly be¬ 
low T° = 0.10. We take m = 4, corresponding to 
q* = (f)0) (see Fig. for model parameters). The 
eigenvector corresponding to the largest eigenvalue of 
M is shown in Table |l] In this case we can simplify 
Eq. (441 by noting that, within the numerical accu- 


Q 

a 

Vaa 

1 

1 

0.489 

2 

1 

-0.701 

3 

1 

0.489 

4 

1 

0.000 

1 

2 

0.000 

2 

2 

0.103 

3 

2 

-0.068 

4 

2 

0.000 

1 

3 

0.000 

2 

3 

-0.017 

3 

3 

0.000 

4 

3 

0.000 

1 

4 

-0.068 

2 

4 

0.103 

3 

4 

0.000 

4 

4 

0.000 


TABLE I: Eigenvector Vaa corresponding to the largest eigen¬ 
value of the pairing kernel M for superconductivity in an axial 
bond-ordered phase. Results are for Jac = 1.31 and rii = 0, 
and the parameters are otherwise as in Fig. [T] 


racy of our calculations, 0.701 « •\/2(0.489), so that 
viiCOs{ka; + f) + vi 3 sm{kx + f) ~ |i; 2 i | cos(A::j,). Fur¬ 
ther, using sin(fca; + f) = cos{kx — f) we obtain. 


Ak(q) = V2A(T)|(5q,ok2i|[cos(A;a;) - cos(fcy)] 
V 32 cos “ 0 + ^22 cos(fcy) 


+^q,2q*W23 COs{ky) 


+S. 


'qi-q* 


Z)14C0S (jta: + ^ -h 1124 COS (fcy)j |. 


where the factor of comes from the definition of 
and we used the equivalence of — q* and {m — l)q*. All 
possible harmonics of q* are present in Ak(q); however, 
the q = 0 component is largest by far and it has a pure 
dx 2 -y 2 symmetry. 

Figurej^b) shows the dependence of the different com¬ 
ponents of Ak(q) as a function of impurity concentration 
for the case Jgc = 1.31. As bond order is reduced by im¬ 
purities, the superconducting components at ±q* and 
2q* make up a progressively smaller fraction of Ak(q). 
When bond order is completely suppressed, these com¬ 
ponents vanish and the system becomes a dirty d-wave 
superconductor. 

The real-space pair amplitudes A^- are more physically 
transparent than Ak(q). Taking nearest neighbor sites i 
and j, 

Ab = (46) 

q k 


where R = {X,Y) = (r^ + rj)/2 and Vij = — vj. We 
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obtain 


A 


j+xj 


V2A{T) 


1 

-^V2l + V32 COS 



(47) 


A 


j+yj 


V2A{T) 


1 

-V21 + V22 cos 



+ ^V23 cos(7rA) 


(48) 


Similarly, the real-space bond order parameter can 
be obtained by inverting Eq. Plots of Ptj and A^ 
with the homogeneous q = 0 component removed are 
given in Figs.j^c) and (d) respectively. These hgures ex¬ 
plicitly show that the spatial modulations of the pairing 
amplitude and bond order are correlated. 


V. DISCUSSION 

The results in Figs. and suggest a way to probe 
possible relationships between charge order and the pseu¬ 
dogap, namely to track the dependence of the pseudo¬ 
gap on zinc doping. We can compare to a number of 
early experiments that explored exactly this, principally 
in YBa 2 Cu 408 , which is often seen as a model under¬ 
doped cuprate because it is stoichiometric. We note, 
however, a well-known and persistent problem that be¬ 
cause the pseudogap appears as a crossover rather than 
a phase transition, the identihcation of the relevant tem¬ 
perature scale(s) depends on the experimental technique, 
and on how the temperature scales are defined. 

Julien et noted that early experiments on 

underdoped YBa 2 Cu 306 +a and on YBa2Cu408 found 
two distinct temperature scales, with dramatically dif¬ 
ferent responses to Zn impurities. The higher scale, 
T* ^ 200-300 K, was seen originally in Knight shift 
measurement^^ that indicated a reduction of available 
spin excitations below T*. Later optical conductiv¬ 
ity measurements showed that there is an accompany¬ 
ing reduction in available charge excitations.^SSEZl Ex¬ 
perimentally, T* was found to be independent of Zn 
concentration.E3EIl The lower temperature scale T' ~ 
150 K wa s observed as a downturn in the NMR relaxation 
rat^Mlll and in the in-plane Hall coefficient.!^ The 

downturn in the Hall coefficient has recently been tied to 
the onset of charge order at Tcp.^^^^Thi s lowe r tempera¬ 
ture is rapidly suppressed by Zn doping.l^Slini _ 

In particular, Zn doping experiment d^^ l ^Q f^ on 
YBa 2 (Cui_iZn 2)408 found that was suppressed from 
TO « 80 K for z = 0 to Tp = 0 K for z ~ 0.03, while 
T' was suppressed much faster,!^ f rom 1 50 K at z = 0 
to 0 K at z « 0.02. Similar resultd^^*^ were found for 
YBa 2 (Cui_^Zn 2)3 06.63- Qualitatively, these are consis¬ 
tent with the suppression of Tbo shown in Fig.[^ To make 
a quantitative comparison, we note that Zn substitutes 
preferentially for Cu sites in the Cu02 planes so that in 


YBa2(Cui_2Znj;)408 the Zn concentration per planar Cu 
is 2z. With this in mind, it is clear that our calculations 
overestimate reduction of superconductivity by disorder, 
relative to experiments. This is a known problem with 
disorder-averaged calculations of Tc in cuprates, which 
neglect spatial inhomogeneity of the order parameter.S3 

Although we have suggested that T' and Tbo may be 
the same temperature scale, we emphasize that a di¬ 
rect comparison between the suppression of T' by Zn in 
cuprates and the suppression of Tbo by impurities in our 
calculations is not straightforward. In particular, Zn im¬ 
purities are known to nucleate magnetic moments locally 
around each impurity site.!^ NMR Ti measurements are 
certainly affected by these moments, and indeed it has 
been suggested that they are sufficient to explain the 
doping dependence of 1/TiT.I^In practice, it may be dif¬ 
ficult to disentangle the contributions of local moments 
and impurity scattering to the suppression of charge or¬ 
der in the cuprates. 

Finally, we remark that the rapid suppression of charge 
order in YBCO by Zn impurities is in cont rast to 
the apparent enhancement of stripe correlation^^^Ml jn 
Zn-doped La 2 -a;Sr 2 ;Cu 04 (this point was also made in 
Ref. [T6|l . We take this as further evidence that the 
physics underlying charge order in YBCO and BSCCO 
is different than that in the La-based cuprates. 


VI. CONCLUSIONS 

We have studied the effects of strong-scattering point- 
like impurities on charge order and superconductivity in 
the cuprate superconductors. Calculations were based on 
a one-band model in which bond order is the analogue 
of charge order in the cuprates. Impurity effects were 
described with a self-consistent t-matrix approximation. 

Our main observation is that d-wave superconductiv¬ 
ity is more robust against impurities than bond order; 
this implies that charge order in the cuprates should be 
more rapidly reduced by Zn substitution than supercon- 
dutivity, even though the onset temperature for charge 
order is higher than Tc. Interestingly, the sensitivity of 
bond order to impurities is not directly connected to the 
symmetry of the order parameter, but occurs because 
charge order arises from only small “hotspot” regions of 
the Fermi surface. 

Experimentally, both the pseudogap and stripe phase 
in cuprate high temperature superconductors are insen¬ 
sitive to Zn doping. This is inconsistent with simple sce¬ 
narios in which charge order contributes directly to the 
pseudogap. 
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Appendix A: Impurities in the bond-ordered phase 

We use the self-consistent t-matrix approximation 
(SCTMA) to obtain an expression for the self energy 
'S{iuJn) due to the impurities. The SCTMA gives the 
disorder-averaged Green’s function and is exact in the 
limit where the impurity concentration rii is small.l^ 
Apart from the complications arising from the charge or¬ 
der, our approach is standard. 

The derivations in this appendix have three parts. In 
Appendix EH the scattering self energy for weak bond 
order is obtained to linear order in Pk(q); this is used 
to obtain the self-consistent equations for Tbo- These are 
solved in Appendix |A 2| to find an approximate analytic 
expression for Tbo. Finally, in Appendix |A3| we find the 
self energy for the case of arbitrarily strong bond order 
with period-TO commensurability. In this case, S(tLi) is 
an m X m matrix. 


ble positions for the /th impurity. To obtain the second 
term in Fig. [^a), we keep only second-order scattering 
contributions in which both impurity lines are from the 


I i _ ) ^ ^ ) I ) ^ 

k” k" k'" k'" k" k" 

• • • 


k" k" + q k'" k'" k" k" + q k" k" + q k'" k'" 

FIG. 8: SCTMA diagrams contributing to the self energy to 
(a) zeroth order and (b) first order in Pk(q). Internal wavevec- 
tors k" and k'" are summed over. Dashed lines represent 
scattering by the impurity potential Vi, and solid lines repre¬ 
sent the Green’s functions G(k", k”; ia;„), G(k"-f q, k"; iLU„), 
etc., as indicated by the momentum labels. 


same impurity. This gives 


1. Linearized Results near Tbo 

In this section, we derive Eq. ( |19| ), along with expres¬ 
sions for the self-energy components and EJi which 
are valid to zeroth and first order in Pk(q) respectively. 
We consider a dilute distribution of Ni pointlike impuri¬ 
ties. Each impurity is assumed to shift the potential on a 
lattice site by Vi, and we will make use of the assumption 
that Tii = NijN <C I, where N is the number of lattice 
sites. The potential energy of electrons interacting with 
the impurities is 

Ni Ni 

F = E ha, = F, ^ E 4.ck'. 

1^1 k,k' 7^1 <7 

_ (Ai) 

where R/ is the position of impurity / and is the 
electron charge density operator on site R/. 

The impurity self energy is obtained by disorder¬ 
averaging over the possible positions R/ of each impu¬ 
rity, and retaining all irreducible diagrams that are first 
order in rii. Figure shows diagrammatic contributions 
to and EJ*. The first term in Fig. |^a) is 

T/. 

^E = niVid^M', (A2) 

^ /=1 


where (...) = ^ X^a/I- • ■] average over all possi¬ 


T /.2 

# E <5z.,EG'(k",k";za;„) 

I,.J=l k" 

^ ^g-*(k-k").Rjg-*(k"-k').Rj^ 

= n,F/5k.k'5° (AS) 

where = N~^ G(k", k"; iujn)- Following this pro¬ 

cedure, the jth order diagram is then ni[V)5)(]-^(5k,k': and 
the sum of diagrams to infinite order is S„Jk.k', where 


E 


n — 


niVi 


(A4) 


The sum of diagrams of the type shown in Fig.|^b) can 
be obtained in similar fashion. There are j — 1 terms at 
jth order in F: each of these terms contains j — 2 factors 
of and one factor of ^^k" G(k"-l-q, k"; ia;„). 

The sum of diagrams is E()(5k'.k+q where 


E3 = 


riiV^G^ 


l + 2Vig°+3V^g°J 

I 

[1 - 


(A5) 


The equations for E„ and E)* are made self-consistent 
by obtaining equations for and G^. These come from 
the equations of motion for the Green’s function. 


G(k,k;fa;„) 

G(k, k-f q;iwn) 

iUJn ^k ^k(Q) 


1 0 

G(k -b q,k;fa;„) 

G(k-b q, k-b q;ia;„) 

^k(Q) ^k+q '^ri 


0 1 


(A6) 
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from which, 
G(k, k;ja;„) 
G(k + q, \s.;iujn) 


1 


(A7) 


_ a(q) + ^ 

{iuJn ^k+q Sn)(^^n ^k ^n) 

(A 8 ) 


to linear order in Pk(q)- 


Equations (A4) and (A7| form a closed set of self- 


consistent equations for E„. Once S„ is known, Eq. (A5) 


and Eq. (A 8 ) can then be solved self-consistently for 

sq =__ (A9) 

with 


= - V 


1 


N 


(iUJn ^n)(^^n ^k+q ^n) 


and 

AT 


a(q) 


A^ ^k ^n)(^^n ^k+q ^n) 


, (AlO) 


• (All) 


Erom Eq. (A9), one can express as 


(A12) 


where P“(q) is defined by Eq. (15). Once E„ and E^ are 
known, we substitute Eq. (A 8 ) into Eq. (11^ to obtain 


p = _ ^k(^k + __ 

N ^ {iWn - Ck+q - E„)(iw„ - Ck - E„) ’ 

(A13) 


k.n 


which is Eq. (19) in the text. In the clean limit (E„ = 
Eq = 0), this reduces to 


Ao/3(q) = 


/3/(ek+q) - /(ek) 


Ck+q — Ek 


(AM) 



FIG. 9: Illustration of the dispersion £(k) near the Fermi 
surface points nested by the bond ordering wavevector q. 
Solid black curves are Fermi surface segments and q* con¬ 
nects hotspot points 1 and 2, at which the joint density of 
states is large. The dispersions in the neighborhood of points 
1 and 2 are e(q) = vpPx + Kpy and e(q + q*) = —vfPx + KPy 
respectively. In this figure, the origin of the coordinate system 
is hotspot 1, not the center of the Brillouin zone. 


and absorb the real part A/x into the chemical potential. 
This approximation is not entirely justified, owing to a 
nearby van Hove singularity in the density of states; how¬ 
ever, we have found that adding a weak linear energy de¬ 
pendence to 7 does not change our answers appreciably. 

For definiteness, we will consider period-m axial or¬ 
der, with Qx = Q* = ‘^T^lrn and qy = 0. We know from 
numerics that only the basis functions 

* 

77 j( = A/2cos(fca;-I- y); r]^ = V2 cos{ky) (A16) 

contribute to Pk(q*), so we restrict our discussion to the 
2x2 subspace in which a,/3 € {1,2}. For axial order, 
hotspot 1 is, in the original coordinate system, at kx = 
2n — q*/2 and ky = n (see, for example. Fig. I^a)), and 
we approximate and 77 ^ by their values at this point. 
We thus obtain 


2. Analytic approximation for Tbo 


ql^V2-, r,l^-V2. (A17) 


We begin with Eq. (19) for the bond ordering kernel 
Eh^(q*) and make a number of simplifications. First, we 
assume that q* nests two Fermi surface hotspots, labelled 
1 and 2 in Fig. that are characterized by anti-parallel 
Fermi velocities vp and by curvatures k. By expanding 
the dispersion around the hotspots we obtain 


Ek = vfPx + i^pI; ek+q* = -vfPx + kpI, (A15) 

where q is the wavevector measured relative to hotspot 
1 (whereas k is relative to the Brillouin zone center). 

Then, we make the approximation that the scattering 
self energy is piecewise constant, so E„ = A/i—f 7 Sgn(a;„) 


Numerically, we find that EJ* is small and, neglecting 
it, we obtain 


Eii(q*) = F 22 (q*) = -F 2 i(q*) = -Fi 2 (q*). (A18) 

with 




q „ Cq+q*)(*W„ Cq) 


, (A19) 


and An = Wn + 7 Sgn(wn). 

According to Eq. ( |I^ , the onset of bond order occurs 
when the largest eigenvalue of F(q*) is equal to 1/Jbo- 
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The eigenvectors of F(q*) are 0 and 2^11 (q*), so Tbo 
satisfies 


l = 2JboFii(q*). (A20) 

The corresponding eigenvector of F(q*) is (1,-1)^, 
which gives the solution -P^(q*) = — P^(q*), 

or 

Fk(q*) = P'(q*)(^i^-ry^), (A21) 


similar to that found numerically. Our goal is now to 
estimate Pii(q*). 

Transforming the summation over q to an integral, 
Eq. (|A19|) becomes 


T 

'^n,ab = I + 

I 

I 

a b 


+ 


A 

' I ' 


+ ... 


dh 


e f db 


FIG. 10: Diagrams used in the SCTMA calculation of the 
self-energy matrix element T.n,ab- Dashed lines represent scat¬ 
tering by the impurity potential Vi, solid lines are the Green’s 
function, and indices a, b, . .. label reduced Brillouin zones for 
the initial and final momentum of the scattered particle. In¬ 
ternal indices are summed over. These diagrams include scat¬ 
tering to all orders from a single impurity. In the SGTMA, 
the Green’s function contains the self energy due to impurity 
scattering, and is calculated self-consistently. 


^- 2 T rr ^ t 

Jo Jo 


{A22) 


{x + iUn - KpI){x - ioJn + Kp^) 

where x = vpPx- The term is a large-energy cutoff, 
and is assumed much bigger than any other energy scale 
in the calculation. 

E valuat ing the integral over x, and substituting into 
Eq. (A20) gives an equation for the bond ordering tem¬ 


perature. 


where K is a reciprocal lattice vector of the original lat¬ 
tice. 

The potential energy of electrons interacting with the 
impurities is given by Eq. (Al). This can be re-written 
as 


E E<v(k-k')4/kv, 


(A25) 


k.k'GBZ' cr 


where dlka is the column vector defined in Eq. (261, 
V(k — k') has matrix elements 


i=4Jbo— 

Jo 


CJn 


n—0 


UJ 


FJpt 


tan 


-1 



(A23) 

where we have dropped a small logarithmic correction 
that vanishes in the limit ttvf ^ Kpy. 

A similar equation holds for the clean limit transition 
temperature provided we replace by a;„. Setting 
these two equations equal to each other, we obtain 


yA r dpy n + I + 7 

tt (n + ^ + + KcpI 

— dPy ^ + I 

F (n+^)^ + koPl 


(A24) 


where = ^/(STrTbo) and ko = k/{2'!tT^^), and A^ = 
vfI2T\jo and Aq = vf/‘2T^^. The cutoffs Ac and Aq 
come from approximating tan“^(a:) ~ (7r/2)0(a: — 1), 
with ©(cc) the Heavyside step function. Performing the 
sum over Matsubara frequencies, we obtain the final re¬ 
sult, Eq. (23). 


3. Commensurate bond order 


Kj,(k-k') = (A26) 

I 

k and k' are now restricted to the reduced Brillouin zone, 
and a, 6 S [1, m]. 

The self energy is obtained from the sum of non cross¬ 
ing irreducible diagrams shown in Fig. |10| These include 
all diagrams to linear order in rii due to scattering from 
the impurity potential: = Ti\ + -I- ..., where S;! 

is jth order in Vi and the subsript n indicates that the 
self energy is evaluated at Matsubara frequency a;„. As 
a result of disorder averaging, all terms depend on k and 
k' only through a term ^k.k' that conserves momentum. 
The first order term is 

Ei,(k,k';zc^„) = (Fc6(k-k')) 

_ W -i(k-|-aq-k'-6q).R/ 

I Ri 
— ‘Ri^ida,bdk.,k.' 

= (A27) 

where (...) = A [...] is the average over all possible 
positions for the Jth impurity. 

Similarly, the irreducible second order term is 


In this section, we derive a set of self-consistent equa¬ 
tions for the response of Pk(q*) to pointlike impurities 
for the case of commensurate period-m bond order. In 
this case, the ordering wavevector q* satisfies mq* = K, = riil^^Jk.k'Sn.ab (A28) 


KabK^' = E - k")Gcd(k")Ed6(k" - k')) 

k" cd 
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where 

1 1 

Gn,ab — ^ ^ ^ GcdiX",'i^n)Sa-b,c-d, (A29) 

N' m 

k"GBZ' c,d—l 


and where it is understood throughout this appendix that 
the Kronecker delta function Sa-b,c-d is satisfied modulo 
TO. We have also explicitly written N = N'm, where N' 
is the number of k-points in a single reduced Brillouin 
zone. For reference, 

Gn,aa Gn'j Gn,a+la Gn J (A30) 

in the limit of small Pk(q*), where and are defined 
in Appendix I A 1 | 

The irreducible third order term is 


^n^ab 




k'GBZ' c,<i=l 


1 


jY'm ^ ^ , w)(5o_c+d-t),/-g 

k"GBZ' /,g=l 

(A31) 


k"GBZ' /,g=l 
m 

Gn,ahGn^hb • 


h=l 


At this point, the pattern is established: the jth order 
term in the series is a matrix product of j factors of the 
TO X TO matrix Qn- We define a t-matrix 

T„ = v,i + v,^gn + K^gnGn + .-. 

= + ^i^riTn 

= [1-V,g^]-^V,, (A32) 

where [.. .]~^ indicates a matrix inverse. Then the self¬ 
energy matrix is 

Sn = UiTn = niVi[l - Vign]~\ (A33) 


To determine the effect of impurity scatt ering on bond 
order, one must simultaneously solve Eq. (A33l for the 
self energy and Eq. (30) for the order parameter. These 
equations are linked by Eq. (32) for the Green’s function. 


Appendix B: Impurities in the Superconducting 
Phase 


1. Tc equations for dirty superconductors 


We briefly review the equations for dirty super¬ 
conductors in the absence of bond order, calculated 
with the SCTM A. Much of this discussion can be found 
elsewher^^ 2 llll and we include it here for completeness. 
In the absence of bond order, Cooper pairs have zero 
center-of-mass momentum and the Hamiltonian is 


ff = 




(a ^ 


Ck 

Ak 


Ckt 

|_Ckt> 


Ak 

-e-k 


J 


(Bl) 


Because of the particle-hole transformation for the spin- 
down component, the impurity potential is 


i/ = ^y 

N ^ 

k,k' 


['^kf E ' 


-i(k-k')-R/ 


' T3 


Ck't 

J 

^-k'^ 


(B2) 

where T 3 is a Pauli matrix in particle-hole space. The 
impurity self-energy given by summing the SCTMA dia¬ 
grams shown in Eig. f^is then 


= n^ [ViTs + V^T3g{iuJn)T3 -I-...] 

= n^V^T3[l-Vig{iu}n)T3]~^ 

= 


where 


GiiUJn) = -^E 

k 


iuJn “h 6 i(; Ak 
Ak ^k 


= gn'TO + gnTl + gnT3, 


(B3) 

(B4) 


(B5) 

(B 6 ) 


and 


lOJr, 


= iuj„, - E“ = iuj„, - 




Ak — Ak -h — Ak — 


(1 - 

n^v^gk 


^B7) 


^k — ^k + — ek 


(1 - _ (y^go) 

n.I4 (1 - V^g?^) 

( 1 - 1^3)2 _ 


2(B8) 

2 (B9) 


Erom the structure of Eq. (B5), one sees that g „ is p ure 
imaginary, while and are real. Equation (B5) ne¬ 
glects terms of order A^, as these are small near Tc- Tc 
is then obtained by solving the linearized equation 


A' 


k — E ffkffk' (C-ky cjt/t) 

k' 


Ak' 

+ Ew ’ 


JscT 

N 


E 

k' ,n 


5k5k' 


(BIO) 


where ^k = 1 for isotropic s-wave superconductors and 
5 k = cos kx — cos ky for d-wave superconductors. In this 
work, numerical results for Tc without bond order are 
generated by solving Eq. (BIO) self-consistently. 

To illustrate the role of each component of the self¬ 
energy, and in particular the anomalous self-energy 
we take the simple case of a band with a constant density 
of states Nq. The components of g{ibjn) are 

1 


6 “ = -vE 


N ^ ujt, 

k 


= -Ar 


de- 


lUJn 


° J < + e- 
= -iTrNoSgn{uJri) 

Gn = -^0 


= 0 . 


(Bll) 

(B12) 
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It then follows that 

, ^ UiTTNoV^ 


and 


Ck = Ck + A^; Afi = 


riiVi 


l + inNoV^r' 


For d-wave superconductors 


= -V 


All 


= 0 , 


(B13) 


(B14) 


(B15) 


so EJj = 0 and At = A^ = A^^. Then, Eq. (BIO) 
becomes 


, _ JscT 

“ ~W~ ^ i 


5k 


(B16) 


Because = (|u;„| + 7 )^, the effect of 7 is to renormalize 
the Matsubara frequencies away from zero, which is qual¬ 
itatively similar to raising the temperature in Eq. (B16). 
Impurities thus impede d-wave superconductivity. 

For isotropic s-wave superconductors Ak = A and 




A 


=-Ao7rT^0(A - |a;„|), (B17) 


where ©(re) is a step function and A is a cutoff that is 
typically of order the Debye frequency. Then, combining 
Eq. (B17), Eq. (B 8 ) and Eq. (B13), we obtain the self- 


consistent equation 


A = A-h 7 pr-r, 


{\uj„\ < A) 


which has the solution 


A 

|w„| 


A 

IWnI 


(B18) 


(B19) 


This result is directly relevant to the Tc equation, which 
in this instance is given by Eq. (BIOI with = 1: 


A = 


No JscT Y, / de-- 

= TT No JscT Y 7^ 


A 




A 


- ttNqJscT Y I 


(B20) 


I^tiI<a 


The last equality follows from Eq. (BI9), and the switch 
of the constraint from |w„| < A to |a;„| < A introduces 
an error ^ 0 ( 7 /A). The key point of this derivation is 
that the the anomalous impurity self-energy S^, which 
renormalizes A, cancels the renormalization of a;„ by , 
so that the Tc equation is the same as in the clean limit. 
In the d-wave case, where = 0, Tc is reduced by im¬ 
purities. 

2. Tc in the bond ordered phase 

In this section, we derive the linearized self-consistent 
equation for the superconducting order parameter Ak(q) 
in the bond ordered phase. From Eq. (38), we have 




A“(aq*) = ^ ^C(c-k,+aq-tCk,t) 


keBZ' e=i 


J S( 

mN' 


E (B21) 


kGBZ' 1=1 

where d^k is the mx m anomalous Green’s function with 
matrix elements 


[‘d'k]a6 (Cka^B_kf,^), 


(B22) 


and where it is understood that £ — a is modulo m. To 
obtain d'k, we solve the equations of motion: 


iujn - Hk(q*) - 
-a! - Sn 


ILOr, 


— Ak — Sn 
H_k(-q*)^-S„ 



G T 


1 0 


T G 


0 1 


(B23) 


to linear order in Ak. To simplify the calculations, we make the approximation that S„ = 0, which is strictly true 
for pure d-wave superconductors. We find in our numerical solutions that the non-d-wave components induced by 
the charge order are typically an order of magnitude smaller than the d-wave components, so that this result remains 
approximately true. Then, we obtain the m x m matrix 


J-k = [zw„l - Hk(q*) - [Ak] [iuJr^l -I- H_k(-q*)^ - S„] 


(B24) 
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Combining this with Eq. (B21), we obtain 




mN' 


n keBZ'e,e'=i p 


r 


(B25) 


This is the result shown in Eq. (|43|). We show in Ap- Thus, 
pendix 


B3 


that 


3. Impurities at Tc in the bond ordered phase 


In the superconducting state, Eq. (A25) gives the po¬ 


tential energy of the impurities in the spin-up block. In 
the spin-down block, we make a particle-hole transforma¬ 
tion and let k —> — k. The particle-hole transformation 
introduces a minus sign, but leaves the form of the poten¬ 
tial otherwise unchanged. Then, combining both spin-up 
electrons and spin-down holes, we obtain 


E 




k.k'GBZ' 


V(k-k') 0 

0 -V(k-k')* 


«'k' 


(B26) 

where is a rank-2m array of particle/hole annihilation 
operators, defined in Eq. (39). 


Because superconductivity modifies and at sec¬ 
ond order in Ak, it is neglected in the linearized equations 
near T^. The equations for and a re thus obtained 
by setting Ak = Sn = 0 in Eq. (B23) for the Green’s 


functions, and then performing the SCTMA sums shown 
in Eig. [T^ Because the particle and hole blocks are de¬ 
coupled in both the Green’s functions and the impurity 
potential, and can be evaluated indep enden tly. 

To linear order in Ak, Sn is given by Eq. (A33l; 
satisfies an equation at similar to Eq. (A33), but with 
kj -t -17 and Qn,ab -t Qn,ab, where 


- TliVi [^1 -f I7^(za^yj)j 


-1 


(B28) 


Because our solutions for Tk(q) involve only cos(A: 2 ; -I- 
qx/‘2) and cos{ky + Qy/^), it follows that H_k(—q) = 
Hk(q) and Hk(q)^ = Hk(q)* (ie. the matrix Hk(q) is 
Hermitian). Then 


.,ab — 


E ^n]cd 

k^BZ' c,rf 

^ ^a—c,b—d 


1 


^ ^ 77k(q) + S: 


mN' ^ 

k^BZ' c,d 
^ ^a—c,b—d 


(B29) 


Substituting this latter form into Eq. (B28), it follows 


that — and satisfy the same self-consistent equa¬ 
tion. We then make the identification 


1 


lb — 


mN' E E ‘^cd(kGa;„)(5a-b,c- 

k^BZ' c^d—1 
- m 

— E E K+H_k(-qr-s 


mN' 

k^BZ' c,d—l 
^ ^a—c,b—d 


^ 1 -1 
cd 


(B27) 


s„ = -s:. 


(B30) 
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